Parity Forbidden Excitations of Sr 2 Cu0 2 Cl 2 Revealed by Optical Third-Harmonic 

Spectroscopy 
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We present the first study of nonlinear optical third harmonic generation (THG) in the 
strongly correlated charge-transfer insulator Sr2Cu02Cl2- For fundamental excitation in the 
near-infrared, the THG spectrum reveals a strongly resonant response for photon energies near 
0.7 eV. Polarization analysis reveals this novel resonance to be only partially accounted for by 
three-photon excitation to the optical charge-transfer exciton, and indicates that an even-parity 
excitation at 2 eV, with a% g symmetry, participates in the third harmonic susceptibility. 



A central theme of condensed matter physics is 
the problem of treating Coulomb interactions on the 
electron- Volt scale to determine such low-energy phe- 
nomena as the Mott metal-insulator transition, superex- 
change, and superconductivity. The two-dimensional 
copper oxide plane has emerged as a model system for 
studying these phenomena. While excitations at the 
high-energy scale have been observed in the the linear op- 
tical spectrum of insulating cuprates, the details of their 
lincshapes, temperature dependence, and even their as- 
signment has been controversial Q]. In this Letter, we 
employ nonlinear optical spectroscopy to determine the 
symmetry and spectrum of both even and odd parity 
states near the charge-transfer (CT) gap of the model 
insulating cuprate, S^CuO^C^. Through polarization 
analysis of the nonlinear optical response, we find ev- 
idence for an even parity excitation at approximately 
2 eV, with a± g symmetry, and we estimate the magni- 
tude of the dipole matrix element that couples this exci- 
tation to the CT exciton. The technique that we describe 
here is generally applicable to correlated insulators, and 
should provide valuable input into models of their elec- 
tronic structure. 

The undoped one-layered oxyhalides of the M2C11O2X2 
(M = Ca, Sr, Ba; X = F, CI, Br) family are of special 
interest for studying the excitation spectrum in cuprates. 
These compounds share the quasi two-dimensional Cu02 
layers which form the structural and electronic basis 
for high-temperature superconductivity. Unlike their 
pure oxide counterparts, however, the intervening alka- 
line earth and halogen elements provide these materials 
with a very stable stoichiometry. Sr2Cu02Cl2 in particu- 
lar can be obtained in single crystal form with very high 
quality H . It is regarded as an almost ideal realization of 
a two-dimensional, spin- 1/2 Heiscnbcrg anti-ferromagnet 
at half filling and serves as the test compound for theory 



in the low-doping regime. 

The most obvious optical feature in the insulating 
cuprates is the CT gap excitation near 2 eV. Cluster cal- 
culations generically predict additional, optically forbid- 
den CT excitations, with a\ g (s) or b\ g (d x 2_ v 2) symme- 
try in A/mmm Q^j, but presently there is insufficient 
experimental guidance for selecting an appropriate min- 
imal model |(| [7], |). Parity forbidden excitations with 
ligand field, or dd character, are also expected; these are 
intraatomic excitations that result directly from strong 
correlations, in which the hole in the ground state (g) 
Cu(3d x 2_ y 2) orbital is transferred to another Cu orbital 
of different symmetry. 

Large energy-shift laser Raman spectroscopy provided 
early evidence for the existence and energy spectrum 
of such excitations, with the assignment of a transition 
d x 2_ y 2 — > d xy at 1.35 eV also identified in optical 
measurements [[Io| [Uj. From a ligand-field analysis of 
related peaks in the resonant X-ray Raman spectrum, 
Kuiper et al. have suggested that the d 3z 2_ r 2 (ai g ) state 
should be located at 1.5 eV, but were unable to observe 
it directly for technical reasons [jl2| . While a peak in the 
Big laser Raman spectrum of some insulting cuprates 
may be consistent with an assignment to this state, this 
peak is only observed in T'-phase materials 

In our experiment, we completely characterize 
the THG third order susceptibility tensor elements, 
X^(— 3lu; lu 7 uj 7 uj), for the Cu-0 {x,y}-plane of 
Sr2Cu02Cl2- In the space group IA/mmm possessed 
by Sr2Cu02Cl2 the nonlinear susceptibility tensor has 

(3) 

two independent elements in the {x, y}-plane, Xxxxx and 
(3) 

Xxxyy, which are in general complex. All tensor elements 

(3) 

related to Xxxyy by permutation of the indices are equal, 
Xyyyy = Xxxxx, and all other tensor elements are zero. 
For a driving optical field E at frequency u>, polarized in 
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FIG. 1: Third harmonic generation in Sr2Cu02Ci2- The left 
axis shows the relative values of I(3u) oc |x**axi:| 2 ' The right 
scale indicates absolute values for Si^CuO^Cb- Dashed line: 
guide to the eye. 



the {x, y}-plane, the nonlinear polarization is given by 

P x (3^) = X^ xx E 3 x +3 X ^ yv E x E 2 y , 

P y (3u) = X x 3 L x E 3 y +3 Xx % y E 2 x E y . (1) 

Throughout, we take x and y to correspond to coordi- 
nates aligned with the crystalline axes of S^CuC^C^, 
with the z-axis normal to the crystal and in the di- 
rection of propagation. It is useful to write the two 
tensor elements in terms of an amplitude and a phase: 
Xxlxx = k xx e za , and Xxxyy = K xy e zlS , with the phase 
difference defined as 8 — (3 — a. 

We first measured the magnitude of the third har- 
monic susceptibility |xxx X x| at room temperature, as a 
function of fundamental photon energy. These and all 
other measurements described here were performed with 
a 250 kHz Tksapphire laser amplifier driving an opti- 
cal parametric amplifier. To account properly for vari- 
ations in pulse duration and mode profile as the laser 
wavelength is tuned, we compared the third harmonic 
intensity I(3w) generated in Sr2Cu02Cl2 with a quartz 
reference. We accounted for the Sr2Cu02Cl2 absorp- 
tion at 3ll> in the limit £ <C A, Ak£ «0, appropriate 
here |l4j. The Sr2Cu02Cl2 samples were cleaved to a 
thickness of £ ~ 100 nm, and oriented crystallographi- 
cally with X-ray diffraction. The quartz reference was a 
thin, £ = 150 /im, c-axis plate, whose absolute nonlinear 
susceptibility, x^xlxx (quartz)— 3 x 10 -14 esu, is known in- 
dependently |H| . The measured relative intensities along 
with the resulting IxWxa;! of Sr2Cu02Cl2 are shown in 
Fig. [j]. The spectrum exhibits a broad « 0.5 eV wide 
resonance at 0.7 eV, varying by more than a factor of 20 
over the explored frequency range. Since the suscepti- 
bility x (3) (-3^; 

lo,lo,uj) involves summation over differ- 
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FIG. 2: (a) Relative phase 8 and (b) relative magnitude p of 
Xiiii and Xxxyy in Sr2Cu02Ci2. The dashed lines are guides 
to the eye. Inset: Example of a retardance measurement: 
P(3cj) (filled circles) and F(3u;) (open circles) recorded at 
lu = 0.89 eV. The dotted lines represent best fits to a sin 2 
function, from which 8 and p are extracted. 



ent intermediate states, the resonance may result from a 
parity allowed, three-photon transition at 3w = 2.1 eV, 
a parity forbidden, two-photon transition at 2lu = 1.4 
eV, or both. All of these possibilities are consistent with 
previously published results Q . 

Our chief experimental results are the spectroscopic 
measurements of the relative amplitude p = n xx /K xy 



("3") 

and the relative phase 8 of Xxxxx, X 



X3) 



shown in 



Figs. |](a) and (b), respectively. These data were ob- 
tained from polarization-sensitive measurements of the 
third harmonic optical intensity along each of the two 
crystalline axes of the sample, 1 x,v (3uj), while varying 
the polarization state of the fundamental laser field. In 
the experimental geometry used here, the incident field 
is initially polarized at 45 degrees to the two crystalline 
axes, then passes through a polarization compensator be- 
fore entering the sample. The natural axes of the com- 
pensator are aligned with those of the crystal, so the 
polarization state upon entering the sample is given by 
E = [xE + yEe tA ]/y/2, where A is the compensator re- 
tardance. Thus, the nonlinear polarizations P x (3u>) and 



P y (3uj) involve controlled mixtures of x 



(3) 



and Xxxyy i 



with A as the control parameter, as may be seen through 
Eq. | 

As the compensator is adjusted, the THG intensities 
F(3w) oc |P x (3w)| 2 and F(3w) oc \P v {3u)\ 2 both oscil- 
late sinusoidally with A. Results from a typical mea- 
surement, for huj = 0.89 eV, are shown in the inset of 
Fig. [| Straightforward manipulation of Eq. |l| tells us 
further that the phase offset in A between these two in- 
tensity curves is 8, and that the minima of l x ' v (3u>, A) are 



3 



offset from zero by K xy (p — 3). From a series of such mea- 
surements at various laser frequencies, we obtain the dis- 
persion of the relative amplitude p(u>) and relative phase 
S(uj), shown in Figs. ||(a) and (b), respectively. These 
data show that p(u>) and S(u>) both undergo clear, step- 
like changes at about 1 eV. Since these are relative quan- 
tities, we can derive the model-independent conclusion 
that as the photon energy crosses fku w 1 eV, the transi- 
tion dipole elements involved in Xxxxx and Xxayv connect 
states with different symmetry. We show below that this 
can only be the result of an excitation with a± g symmetry 
at 2 eV. 

The tensor elements of the nonlinear susceptibility are 
given by 

N 8 

(2) 

where the intrinsic permutation operator Vi applies to 
the last three tensor indices {jkl} [Gil. The eight terms 
Ky 1 ' correspond to different summations over intermedi- 
ate states. The dominant term in our experiment is KS 1 ' , 
which may be written 

£■(1) _ f 1 gi> 1 l J 'i 1 L> 2 Vi 2 i> 3 f J i 3 g ,„s 

where {^1,^2,^3} are the intermediate states, p\, lV2 — 
{vi\p, l \v2) are dipole matrix elements along the direction 
f{, and each Q, Vl g — u) Ueg — V)v t g is a complex frequency 
associated with a transition from the ground state g to an 
intermediate state Vj . Individual terms in Eq. [| may be 
one-, two-, or three-photon resonant, or it may contain 
double or triple resonances, depending on the location 
and symmetry of the states vg. If the product of matrix 
elements in Eq. ^| is to be nonzero, symmetry requires 
both that v\ and z/3 transform as e u in A/mmm, and that 
V2 transforms as a\ g (z 2 ), a,2 g (x 3 y — xy 3 ), b\ g (x 2 — y 2 ) 
or &2g (xy). Moreover, the summation over all possible 
permutations in Eq. |^ leads to additional cancellations, 
so that those terms in which V2 transforms as a,2 g or 629 
do not contribute to the THG signal. Consequently, the 
even-parity states which contribute to the THG signal 
are exclusively those with a\ g or b\ g symmetry. 

To evaluate the role of these intermediate states in the 
measured third harmonic spectrum, we compare our re- 
sults from Sr2Cu02Cl2 to the results of a phenomeno- 
logical, independent-level model, commonly used in non- 
linear optics |H| . We assume that the excited states are 
independent levels, each with its own energy, dephasing 
rate, and symmetry. To obtain a realistic estimate of 
the matrix elements, excited state energies, and decay 
rates of the odd-parity states, observed in conventional 
linear optics, we performed a fit of four lorentzian func- 
tions to the experimentally determined absorption spec- 
trum [|l7j G^8|. By attributing each of these lorentzians 



to a different state rj n , where n ranges from one to four, 
we obtain from the fit the energies E n , dephasing rates 
7„, and dipole matrix elements p x gn connecting each state 
r\ n to the ground state. These parameters are listed in 
Table |. We emphasize that the parameters used here 
are chosen merely to yield a good description of the lin- 
ear susceptibility, and subsequently remain fixed for the 
nonlinear susceptibility calculation using Eq. ^|. 

The results of this model calculation are shown as solid 
lines in Fig. 0. Panel (a) shows that the absolute mag- 
nitude of Xxxxx exhibits a peak near 0.7 eV, roughly one 
third of the energy of the first odd-parity excited state 
771 of the model, as expected from a three-photon res- 
onance. Panels (b) and (c) show the that the relative 
amplitude and phase are completely flat as a function of 
frequency, with the value of the amplitude ratio p fixed at 
three, while the phase difference S is fixed at zero. This 
is the behavior expected of a spherically symmetric non- 
linear hyperpolarizability [|l9) . It should be noted that in 
this case, the relative quantities p and 5 are entirely in- 
sensitive to the details of the model parameters, because 
their values are determined by symmetry. We find that in 
general, if the ground state possesses bi g symmetry, then 
any set of excited states containing only b\ g and e u sym- 
metries will display this overall, spherically symmetric 
hyperpolarizability, over the entire frequency spectrum. 
The model, however, does give the right order of magni- 
tude for |Xrara|- In the case of Sr 2 Cu02Cl2, then, the 
magnitude of the nonlinear susceptibility is dominated 
by a three-photon resonance to CT excited states, and is 
only weakly enhanced by two-photon resonances to even- 
parity states [^0] . 

While the overall magnitude and shape of the reso- 
(3) 

nance in Xxxxx is explained well with only the optically 
allowed CT excited states, the structure that we ob- 
serve in S and K xx /n xy implies the involvement of ex- 
citations with different symmetries. As we have seen, 
the only other states that may contribute to the THG 
signal are those with a\ g symmetry. Moreover, when 
these states are included in Eq. ||, the transitions to 
the ai g and b\ g intermediate states will interfere. Since 

(x\x\a lg ) = (y\y\a lg ), and (x\x\b lg ) = -{y\y\b lg ), this 

(3) 

interference will exhibit different behavior in Xxxxx and 
(3) 

Xxxyy, to produce exactly the deviations from spherical 



TABLE I: Energies E n , dephasing rates 7„, and matrix ele- 
ments fig n , used in the model calculations described in 
the text. 
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FIG. 3: Comparison of two model calculations for x ■ ( a ) 
magnitude, (b) relative amplitude a, and (c) relative phase 8. 
The inset shows the level scheme assumed in the model, with 
\Xi) the two-photon states of interest. In all panels, solid 
lines indicate that only the virtual excitation to |Xi) = |6i s ) 
is included. Dashed lines indicate the addition of a resonant 
state at 2.0 eV, \X 2 ) = |o ls ). 



symmetry that we observe. Thus, we extend our model 
calculation by adding an excitation with a,x g symmetry 
at 2.0 eV. The additional matrix elements fi^ n coupling 
this excited state to the four odd-parity excited states 
are given in Table [|. Remarkably, the results shown as 
dashed lines in Figs. || (b,c) are in close agreement with 
the experiment, both in sign and in magnitude. There 
is considerable uncertainty in the parameters obtained 
from this simplified model, so the error on our estimate 
of the energy of the a\ g excitation may be as large as 
0.4 eV. However, even with this error, the midinfrared 
range is clearly excluded |ll| , and the symmetry of the 
state is unambiguously determined through this analysis. 

Finally, we turn to the assignment of this excitation. 
Symmetry-allowed alternatives include a CT exciton, a 
Cu dd excitation, and coupled modes, such as an exciton- 
phonon (i.e., ai g € e u ® e u ). It is likely, however, 
that all of these excitations are active in this range, 
and are strongly mixed. This proposal is supported by 
the anomalously broad Urbach tail observed at low tem- 
peratures in Sr2Cu02Cl2, which has been explained re- 
cently by a model that involves strong coupling between 
the charge-transfer gap excitation and another, optically 
forbidden, excitation nearby in energy fll8|| . Electronic 
structure calculations place both the Cu dd and the a\ g 
CT excitons near 2 eV §, I H, and both modes would 
be expected to couple strongly to the e u CT exciton via 
phonons. As our theoretical understanding of these ex- 
citations grows, our estimate of the matrix elements pf &n 
may provide a further test of this scenario. 

In summary, we have determined the spectrum of 
the full THG nonlinear optical susceptibility tensor for 
Sr 2 Cu02Cl2, using both absolute measurements against 



a standard reference and relative measurements based on 
polarization analysis. The dominant feature in the ab- 
solute measurement is a peak at Hco ~ 0.7 eV, which 
can be explained as a three-photon resonance to the 
CT gap. The relative measurements show that an ad- 
ditional, optically forbidden state near 2 eV participates 
in the THG spectrum, and that this state possesses a\ g 
symmetry. Both the symmetry and the energy of this 
state are consistent with a mixed mode involving both a 
d x 2_ y 2 — > d^ z 2_ r 2 transition and a coupled CT-exciton- 
phonon mode. These measurements were performed with 
fundamental frequencies well away from resonance, in 
contrast to previous Raman work, and the relative mea- 
surements have enabled the estimation of the optical ma- 
trix elements coupling the even-parity state to the CT 
gap excitation. In general, both Raman and THG spec- 
troscopy are different cuts through the three-dimensional 
frequency space that x^ ma Y probe. Further devel- 
opment of nonlinear optical spectroscopy will provide a 
more complete determination of these and other optical 
excitations. 
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